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Human  Body  Regional  Convective  Heat  Transfer  Determination 
Using  Sublimating  Naphthalene  Disks 


Stephen  KW.  Chang  and  Richard  R.  Gonzalez,  U.  S.  Army  Research  Institute  of 
Environmental  Medicine,  Natick,  MA.  01760-5007 


Abstract 

The  regional  convective  heat  transfer  coefficients  (hc)  on  the  human  body 
were  determined  using  sublimating  naphthalene  disks.  Circular  naphthalene  disks 
were  affixed  to  various  body  segments  of  a  stationary,  life-size  manikin,  under 
constant  temperature  and  wind  speed  in  an  environmental  chamber.  The  amount 
of  naphthalene  weight  loss  through  sublimation  was  translated  to  hc  using  the 
Chilton-Colbum  j-factor  analogy  between  heat  and  mass  transfer.  The  regional 
convective  heat  transfer  coefficients  can  be  determined  by  using  strictly  the  heat- 
mass  transfer  analogy,  excluding  any  supplementary  technique  of  cylindrical  body 
segment  approximations  or  other  shape  extrapolations.  The  logarithmic  mean 
density  factor  for  naphthalene  sublimating  in  air  (P^n)  was  also  determined. 
Pam.*  for  the  naphthalene-air  sublimation  environment  is  only  one  third  of  the 
water  vapor-air  diffusion  environment  (Pam)-  P AM,n  is  an  essential  factor  for 
extracting  the  correct  hc  value  from  the  naphthalene  mass  transfer  data. 
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Nomenclature 


Cp  =  specific  heat  (J»kg'1»K1) 

Dv  =  mass  diffusivity  (m2/s) 

hc  =  convective  heat  transfer  coefficient  (W^m  ^K1) 

1^  =  naphthalene  mass  transfer  coefficient  (m/s) 

=  Chilton-Colbum  analogy  j-factor  for  heat  transfer 
jm  =  Chilton-Colburn  analogy  j-factor  for  mass  transfer 
m  =  naphthalene  sublimation  weight  loss  rate  (kg'm^s1) 

Pa  =  naphthalene  vapor  pressure  in  air  (assumed  =  0) 

Pam  =  logarithmic  mean  density  factor  for  water  vapor  diffusion  in  air  (ND) 
pAM,n=  logarithmic  mean  density  factor  for  naphthalene  sublimating  in  air  (ND) 
Pr  =  Prandtl’s  number  (ND) 

Pg  =  naphthalene  surface  vapor  pressure  (Torr  or  mmHg) 

R  =  naphthalene  gas  constant  (0.487  mmHg*m3«kg'1»K1) 

Sc  =  Schmidt’s  number  (ND) 

Ta  =  ambient  temperature  (°K) 
u  =  air  velocity  (m) 
k  =  thermal  conductivity  (W'm^K1) 
p  =  viscosity  (kg»m'«s'1) 

p  =  density  (kg/m3)  _ _ 

i  A cannot on  For  _ 

I  K T 1 3  I 

]  n-ri':  T.V.  □ 

i  r:.  .iro-K.rrd  □ 


Introduction 

The  Chilton-Colbum  j-factor  analogy  between  heat  and  mass  transfer  is  a  well 
established  theory  [ASHRAE  Fundamentals,  1985].  Since  convective  heat  transfer 
is  not  a  directly  measurable  quantity,  degree  and  extend  of  mass  transfer  of  a 
sublimating  substance  (e.g.  naphthalene,  para-dichlorobenzene)  enables  the  accurate 
and  direct  determination  of  the  convective  heat  transfer  coefficient  hc  [Sogin,  1958; 
Kreith  et  al,  1959;  Neal,  1975;  Sparrow  and  Tien,  1977].  The  convective  heat 
transfer  can  then  be  computed  from  hc,  incorporating  pi  the  surface-to-ambient- 
air  temperature  difference.  Extensive  work  in  applying  the  heat-mass  transfer 
analog)'  in  measuring  hc  for  the  human  body  was  reported  by  Nishi  and  Gagge 
[1970a,  1970b].  They  attached  naphthalene  spheres  to  different  body  segments  on 
human  subjects.  The  convective  transfer  of  the  sphere  was  then  related  to  the 
corresponding  body  segment  by  approximating  the  body  segment  as  a  cylinder  and 
using  known  convective  relationship  between  ball  and  cylinder.  In  the  present 
study,  we  used  naphthalene  disks,  and  applied  only  the  heat-mass  transfer  analog)’ 
without  any  body  segment  shape  supplementation.  Since  the  naphthalene  disks 
sat  directly  on  body  segments,  the  regional  hc  was  measured  directly. 
Extrapolation  between  different  shapes  was  unnecessary.  In  applying  the  heat- 
mass  analogy,  Nishi  and  Gagge  bypassed  the  Chilton-Colbu!  , -factors  by 
substituting  Hilpert’s  numerical  data  for  cylinders  directly  into  the  heat-mass 
analogy  equation.  We  followed  directly  and  strictly  the  Chilton-Colbum  j-factor 
equations  in  developing  the  present  theoretical  basis. 
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Theory  of  Heat  -  Mass  Transfer  Analogy 

Mass  transfer  of  naphthalene  sublimation  h^  can  be  expressed  as  [Nishi  and 
Gagge  1970a] 

I'm  =  R  •  Ta  •  A  '  <Ps  -  Pa>  111 


Assuming  the  heat  of  sublimation  is  negligible,  Pg  may  be  considered  as  equal  to 
the  saturated  vapor  pressure  at  Ta  [Sherwood  and  Trass  1960], 


log  Pc  =  11.55  -  3765/T 
10  s  a 

or, 

P  _  (11.55— 3765/T  ) 

»  “  10  a 


The  Chilton-Colbum  analogy  j-factor  [ASHRAE  Fundamentals  1985] 


for  heat  transfer  = - (Pr)^,  with  Pr  = 

11  P*Cp»U 


cp*P 


(3) 


for  mass  transfer  jm  = 


PAM,n«hm  2/n  p 

-  (Sc)^,  with  Sc  = 


u 


p.D, 


{4) 


4 


Now,  equating  the  heat  and  mass  transfer  j-factors, 


h„ 


P*cp*u 


m 


u 


(Sc) 


2/3 


c  (  Sc  )2/3 
Pr  ’ 


(5) 


As  an  example,  consider  the  physical  properties  of  naphthalene  vapor  transfer  at 
30°C  (typical  values): 

k  =  0.02636  W-m'-K1  Dy  =  6.330*  106  m2/s 

Cp  =  1006  J’kg^K'1  p  =  1.165  kg/m3 

p  =  1.869*  10'5  kg.m'.s4 

Pr  =  — P  ^  =  0.7133  Sc  = - =  2.535 

K  P*Dv 

pAM.n  =  0.319  (for  0.18  m/s  <  u  <  1.8  m/s) 


Equation  {5}  then  becomes  h_  =  870.6  •  hw 

v  m 


(6) 
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Experiment 

Circular  naphthalene  disks  were  attached  to  the  surface  of  various  body 
segments  on  a  stationary,  life-size  (body  surface  area  1.68  m2)  copper  manikin  in 
an  environmental  chamber.  Figure  1  gives  a  schematic  diagram  of  the  manikin 
and  shows  the  positions  of  the  naphthalene  disks  on  the  upper  arm,  lower  arm, 
thigh,  calf  and  chest.  A  manikin,  rather  than  human  subjects,  was  used  as  the 
appropriate  model  to  avoid  the  problem  of  perspiration  and  thereby  eliminate  any 
evaporative  contribution  usually  involved  with  human  subjects. 

The  environmental  chamber  was  set  at  an  ambient  temperature  (Tfl)  of  30°C, 
with  dew  point  temperature  at  5°C.  The  regional  temperature  and  air  velocity 
were  measured  at  five  sites:  upper  arm,  lower  arm,  thigh,  calf,  and  chest.  Figure 
1  gives  a  schematic  representation  of  the  locations  of  naphthalene  disks. 
Omnidirectional  thermal  anemometers  (range  0  -  3  m/s)  and  thermistor 
temperature  probes  were  placed  approximately  2  cm  above  and  2  cm  away  from 
the  naphthalene  disk,  in  such  a  way  as  not  to  disturb  the  impinging  airflow  to  the 
disk.  Airflow  in  the  chamber  was  directed  normally  at  the  disks.  The  duration 
of  each  experiment  was  55  minutes. 

Scintillation  grade  naphthalene  with  a  melting  point  of  approximately  80°C 
was  melted  and  then  poured  into  aluminum  disk  casting  cassettes.  The  casting 
quickly  hardens  at  room  temperature,  after  which  a  nonadhering  cover  plate  was 
removed  to  reveal  a  smooth  naphthalene  surface.  Immediately  after  casting,  the 
disks  were  stored  individually  in  airtight  containers.  All  naphthalene  disks  were 
allowed  to  equilibrate  in  the  chamber  at  30°C  for  24  hours  before  using.  For  the 


6 


study,  the  disks  were  first  mounted  into  vinyl  retainers,  which  were  in  turn 
fastened  to  the  manikin  using  Velcro  straps.  Disks  were  weighed  immediately 
before  and  after  an  experiment  on  a  precision  balance  sensitive  to  ±0.01  mg. 

To  determine  the  logarithmic  mean  density  factor  of  naphthalene  sublimating 
in  air,  P^n  in  equations  {4}  and  {5},  a  mass  transfer  experiment  similar  to  that 
performed  by  Nishi  and  Gagge  [1970a]  was  employed.  A  1.21  m  long  magnesium 
bar  was  rotated  at  its  center  by  a  stepping  electric  motor  at  10  constant  speeds 
ranging  from  3  to  30  rpm,  in  a  still-air  chamber.  One  flat  naphthalene  disk  was 
attached  to  each  end  of  the  magnesium  bar,  in  such  a  way  that  as  the  bar  rotates 
the  naphthalene  disks  faced  normally  into  the  air.  The  rotation  resulted  in 
translational  speeds  from  0.183  to  1.83  m/s  for  the  disks.  Since  the  chamber  had 
still-air,  the  translational  speeds  were  also  exactly  the  air  velocities  that  the 
naphthalene  disks  were  facing.  Duration  of  each  experiment  wlS  again  55 
minutes.  From  the  weight  loss  data  of  the  rotating  naphthalene  disks,  hc  was 
computed  using  equation  {5},  but  without  the  P^n  factor,  yielding  in  effect  hp/P^  n. 
When  this  rotating  disk  hjTP^  n  was  compared  to  true  hc  data  of  comparable  wind 
speed  (from  Nishi  and  Gagge  [1970a]),  a  value  for  P^n  was  then  obtained. 

RESULTS 

Table  1  shows  the  regional  air  velocities  and  local  temperatures  surrounding 
the  stationary  manikin.  This  level  of  regional  wind  speed  was  maintained 
throughout  the  experiment.  The  local  temperatures  shown  were  temperature  data 
averaged  over  the  entire  duration  of  the  study. 
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The  naphthalene  disk  cassettes  were  appropriately  curved  to  conform  to 
surface  curvatures  of  the  upper  arm,  lower  arm,  thigh,  and  calf  of  the  manikin. 
Casting  cassettes  used  for  the  chest  were  not  curved.  The  curved  disks  expose 
an  elliptic  rather  than  circular  surface.  The  elliptic  surface  area  of  each  cassette 
disk  was  measured,  shown  in  Table  2,  and  properly  quantified  (during  evaluation 
of  the  m  term  in  Equation  {1}). 

In  Table  3,  hc  was  calculated  using  equation  {5},  without  the  P^  n  factor,  and 
subsequently  compared  to  data  calculated  using  Nishi  and  Gagge’s  predictive 
equation  [1970a]. 

With  the  correction  factor  P^n,  the  convective  transfer  coefficient  h£  was 
computed  and  the  results  presented  in  Table  4.  The  h£  shown  represent  regional 
values  at  the  specific  sites.  For  each  study  run,  hc  was  averaged  over  the  55 
minute  experiment  period.  The  results  in  Table  4  are  the  average  of  15  runs. 

DISCUSSION 

In  order  to  facilitate  the  rate  of  naphthalene  sublimation,  the  environmental 
chamber  temperature  needed  to  be  maintained  as  high  as  tolerable.  In  this  study, 
30°C  was  chosen.  At  this  temperature,  most  human  sulgects  will,  to  a  greater  or 
lesser  extend,  perspire.  Once  perspiration  occurs,  it  will  invariably  contaminate  the 
naphthalene,  and  render  any  naphthalene  weight  measurement  useless.  To  avoid 
the  problem  of  perspiration  contamination,  a  life-size,  stationary  copper  manikin 
was  used.  Since  we  were  interested  only  in  the  convective  heat  loss,  a  manikin 
also  eliminated  any  confounding  evaporative  contribution  u*ually  involved  with 
human  subjects. 
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The  logarithmic  mean  density  factor  of  naphthalene  sublimating  in  air  n 
is  a  necessary  correction  factor  which  results  from  difffusion  of  air  toward  the 
naphthalene  surface  due  to  an  air  partial  pressure  gradient  [ASHRAE 
Fundamentals,  1985;  Fobelets  and  Gagge,  1988],  Conventionlly,  for  the  water 
vapor  -  air  interface,  P^v,  is  approximately  equal  to  unity  and  thus  often  ignored. 
For  the  naphthalene  -  air  interface,  our  results  showed  that  this  value  is  very 
different  from  unity.  To  determine  P^  n,  hc  was  first  computed  using  equation  (5) 
sans  the  P^n  factor,  then  compared  to  known  hc  data  of  comparable  air  velocity 
environment.  In  Table  3,  hc  /  P^n  in  column  [1]  was  computed  using  the 
naphthalene  mass  loss  data  from  the  rotating  arm  mass  transfer  experiment.  For 
comparable  data,  a  Nishi  and  Gagge  [1970a]  predictive  equation  was  employed. 
From  analysis  of  their  extensive  naphthalene-h£  study  data,  Nishi  and  Gagge 
obtained  the  following  power  function  regression  equation: 


hc  =  8.60.U0531 


(7) 


where  u  denotes  air  velocity  in  m/s,  and  0.2  <  u  <  2.0  m/s.  Equation  {7}  enables 
prediction  of  hc  from  air  velocity  measurements.  Column  [2]  in  Table  3  resulted 
from  operating  Equation  {7}  on  the  rotating  arm  air  velocity  data.  Comparison 
of  the  two  columns  of  data  yielded  the  P^  n  values  in  column  [3],  for  the  ten 
rotational  speeds.  P^.n  is  apparently  a  function  of  air  velocity.  Following  Nishi 
and  Gagge’s  example  using  power  function  equation,  Equation  {8}  was  obtained 
from  regression  analysis  (^=0.94). 


AM,n 


=  0.326  •  u 


0.126 


0.18  m/s  <  u  <  1.8  m/s 


{8} 
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However,  because  air  velocity  does  not  appear  as  a  variable  in  the  heat-mass 
transfer  analogy  Equation  {5},  Pam.ii  as  a  function  of  air  velocity  cannot  be 
substituted  in  directly.  An  average  value  PxM.n  was  used  instead  in  Equation  {5}. 
From  the  results  in  Table  3,  an  average  Pam.ii  of  0.319  ±  0.0295  was  obtained. 
Obviously,  the  application  of  this  average  value  is  limited  by  0.183  <  u  <  1.83  m/s. 
Notice  in  Table  3,  for  the  regional  air  velocities  encountered  on  the  manikin,  0.4  - 
0.7  m/s  (from  Table  1),  the  average  P^^  value  is  very  close  to  the  actual  computed 
Pam.ii-  For  example,  at  0.548  m/s  (9  rpm),  the  actual  computed  Pam.ii  was  0.311. 
Although  not  presented  explicitly  as  logarithmic  mean  density  factor  P^  n,  Nishi 
and  Gagge  [1970b!  also  gave  an  indication  that  a  new  heat  mass  transfer 
correction  factor  of  0.37  (1/2.7)  was  essential  when  operating  at  the  naphthalene  - 
air  interface  (see  footnote  of  Table  1  in  Nishi  and  Gagge  [1970b]). 

In  Table  3,  the  range  of  P^.n  fo  the  naphthalene  -  air  interface  is  only 
approximately  one  third  of  the  often  neglected  Pam  (*=1-0)  of  the  water  vapor  -  air 
interface.  Therefore,  unknowingly  assuming  a  unity  value  P^n  for  the 
naphthalene  -  air  environment  will  result  in  substantial  overestimation  of  the 
derived  convective  heat  transfer  coeffiecient  hc.  This  type  of  overestimation  was 
reported  by  Sparrow  and  Tien  [1977],  and  evident  in  an  initial  report  on  h£ 
determination  using  an  articulated  manikin  by  Chang  et  al  [1988]. 

Incorporating  the  mean  density  factor  Pam.ii  (average  value  0.319),  regional  hc 
for  the  stationary  manikin  was  computed  and  is  presented  in  Table  4.  The  range 
of  hc  was  between  5  -  8  W/(m2»K).  Since  the  naphthalene  disks  conformed  to  the 
body  segment  curvatures  and  sat  directly  over  the  specific  body  sites,  local  hc  over 
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individual  body  sites  were  measured,  rather  than  some  average  hc  values  for  entire 
body  segments. 


CONCLUSION 

Using  the  sublimating  naphthalene  disks  technique,  the  regional  convective 
heat  transfer  coefficients  (h£)  on  the  human  body  were  determined  on  a  stationary, 
life-size  manikin.  Based  on  the  Chilton-Colbum  j-factor  heat-mass  analogy,  the 
sublimating  naphthalene  disks  becomes  a  straightforward  technique  for  determining 
the  convective  heat  transfer  coefficient.  This  analogy  as  shown  in  this  report 
eliminates  the  many  very  elaborate  ex'  rapolations  employed  by  other  techniques. 
The  need  to  represent  the  human  body  as  an  assembly  of  cylindrical  segments  of 
various  sizes  is  thereby  unnecessary.  The  necessity  of  determining  characteristic 
diameters  for  the  representative  cylindrical  segments,  whose  dimensions  are  critical 
in  Nishi  and  Gagge’s  formulation,  is  eliminated.  The  extensive  analogies  applied 
using  the  relationship  between  the  various  nondimensional  quantities  of  heat 
transfer,  e.g.  Nusselt,  Prandtl,  Reynolds,  Schmidt,  Sherwood,  and  Stanton  numbers, 
are  also  reduced. 

In  this  report,  the  logarithmic  mean  density  factor  for  naphthalene 
sublimating  in  air  (Pam.ii)  was  also  determined.  Pam.ii  for  the  naphthalene  -air 
interface  is  only  approximately  one  third  of  the  often  neglected  Pam  (=1.0)  of  the 
water  vapor  -  air  interface.  P^n  is  apparently  a  function  of  air  velocity.  An 
regression  equation:  Pam.h  =  0.326  •  u0,126  was  calculated  (r2=0.94)  for  0.18  m/s  < 
u  <  1.8  m/s.  Because  the  Chilton-Colbum  heat-mass  transfer  analogy  equation 
does  not  operate  on  air  velocity  as  a  variable,  an  average  value  Pam.ii  =  0.319  was 
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used  in  computing  hc.  P^  n  is  an  essential  property  for  extracting  the  correct  hc 
value  from  the  naphthalene  mass  transfer  measurement  data.  An  overestimation 
of  the  convective  heat  transfer  coefficient  hc  results  without  suitable  assessment 
of  this  factor. 
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Table  1 


[S.KW.  Chang  and  R.R.  Gonzalez] 


naphthalene 

disk 

site 

regional 
air  velocity 

(m/s) 

temperature 

(°C) 

Upper  Arm 

0.61  ±  0.02 

29.97  ±  0.11 

Lower  Arm 

0.71  ±  0.01 

29.91  ±  0.12 

Thigh 

0.42  ±  0.01 

30.00  ±  0.13 

Calf 

i 

0.69  ±  0.02 

29.95  ±  0.11 

Chest 

0.49  ±  0.01 

30.00  ±  0.13 

Table  2 


[S.KW.  Chang  and  R.R.  Gonzalez] 


site 

disk 

surface  area 
(cm2) 

Upper  Arm 

19.23 

Lower  Arm 

19.32 

Thigh 

18.85 

Calf 

19.32 

Chest 

18.09 
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Table  3 


[S.KW.  Chang  and  R.R.  Gonzalez] 


[i] 

[2] 

[3] 

air  velocity  II  from  Eqn.  5 

U  (m/s)  ||  hC/PAM,n 

II 

from  Nishi  &  Gagge 

ratio  [2]/[l] 

rpm 

hc=8.60.u0531 

PaM.I1 

3 

0.183 

13.76 

3.46 

0.251 

6 

0.365 

17.64 

5.07 

0.287 

9 

0.548 

20.16 

6.26 

0.311 

12 

0.730 

22.75 

7.28 

0.320 

15 

0.913 

24.89 

8.18 

0.329 

18 

1.10 

27.02 

9.09 

0.336 

21 

1.28 

i 

29.01 

9.59 

0.331 

24 

1.46 

30.83 

10.51 

0.341 

27 

1.64 

32.79 

11.18 

0.341 

30 

1.83 

34.48 

11.85 

0.344 

average  0.319  ±  0.0295 

0.126 


regression  equation:  Pam  ,, 


=  0.326  •  u 


Table  4 


[S.KW.  Chang  and  R.R.  Gonzalez] 


naphthalene 

disk 

site 

local  hc 

W/(m2»K) 

Upper  Arm 

7.49  ±  0.133 

Lower  Arm 

7.23  ±  0.122 

Thigh 

6.01  ±  0.110 

Calf 

7.78  ±  0.165 

Chest 

5.30  ±  0.120 
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Figure 

Figure 

Table  1 

Table  2 
Table  3 

Table  4 


&  Table  Captions 

Schematic  diagram  showing  locations  of  the  naphthalene  disks  on  the 
manikin 

Regional  air  velocities  and  local  temperature  surrounding  the  stationary 
manikin 

Naphthalene  disk  surface  areas 

Determination  of  P^n  from  the  rotating  arm  mass  transfer  experiment 
radius  =  58.1  cm  (from  center  of  bar  to  center  of  mounted  disk) 
hc  is  in  W'm^K'1 

Regional  hc  on  the  stationary  manikin 
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